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In this paper we consider in detail the properties of dynamical heterogeneity in lattice glass 
models (LGMs). LGMs are lattice models whose dynamical rules are based on thermodynamic, as 
opposed to purely kinetic, considerations. We devise a LGM that is not prone to crystallization and 
displays properties of a fragile glass-forming liquid. Particle motion in this model tends to be locally 
anisotropic on intermediate time scales even though the rules governing the model are isotropic. 
The model demonstrates violations of the Stokes-Einstein relation and the growth of various length 
scales associated with dynamical heterogeneity. We discuss future avenues of research comparing 
the predictions of lattice glass models and kinetically constrained models to atomistic systems. 



I. INTRODUCTION 



The cause of the dramatic slowing of dynamics close 
to the empirically defined glass transition is a subject of 
great continued interest and debate 0, • Different theo- 
retical proposals have been put forward aimed at describ- 
ing some or all of the phenomena commonly observed in 
exp eriments and computer simulations [3|, |J, [5|, [6|, [Tj, |8|, [9|, 
[Mill El El- While these proposals are often based on 
completely divergent viewpoints, many of them are able 
to rationalize the same observed behaviors. This fact 
stems from the somewhat limited amount of information 
available from experiments and simulations. Since the 
growth of relaxation times in glassy systems is precipi- 
tous, it is very difficult, and in some cases impossible, to 
distinguish models solely on the basis of different predic- 
tions of gross temperature dependent relaxation behav- 
ior. In addition, computer simulations, which are often 
more detailed than experiments, are limited by the range 
of times scales and sizes of systems that can be studied. 
These difficulties have hindered the search for a consen- 
sus on the microscopic underpinnings of vitrification. 

Despite the continued debate that revolves around the 
theoretical description of supercooled liquids and glasses, 
little argument exists regarding the importance of dy- 
namical heterogeneity as a key feature of glassy behavior 
, [IE [I?], O, [HI • Dynamical heterogeneity refers to 
the fact that as a liquid is supercooled, dynamics become 
starkly spatially heterogeneous, requiring the cooperative 
motion of groups of particles for relaxation to occur. Dy- 
namical heterogeneous motion manifests in several ways, 
and leads to violations of the Stokes-Einstein relation 
[20L I2H I22I l23l . [2^ | , cooperative hopping motion reflected 
in nearly exponential tails in particle displacement func- 
tions [E HE HE H3, Hl| , and growing length scales such 
as those associated with the recovery of Fickian diffusion 
, growing multi-point correlation functions 
JllH (HI m, H3, S IE EE El. Indeed, the rela- 



tively recent explication of the phenomena of dynamical 
heterogeneity has dramatically shifted the focus of the 
field and has placed new constraints on the necessary in- 



gredients for a successful theory of glass formation. 

Given the similarity of some aspects of dynamical het- 
erogeneity to critical fluctuations in standard critical phe- 
nomena, it is natural to investigate two and three dimen- 
sional simplified coarse-grained models that encode the 
crucial features of this heterogeneity. Currently, the most 
investigated class of coarse-grained models are the "ki- 
netically constrained models" (kcms) H, [EEEEEE3- 

KCMs are spin or lattice models that generate slow, 
glassy relaxation via constraints on the dynamical moves 
that are allowed. The slowing down of the dynamics is 
caused by rarefactions of facilitating regions, also called 
defects. Importantly, although the dynamics is com- 
plex the thermodynamics is trivial since the dynamical 
rules are such that all configurations are equally likely. 
The philosophy of this viewpoint is that thermodynamic 
quantities, such as the configurational entropy, are not 
the fundamental underlying cause of the growing time 
scales in supercooled liquids. It has been argued that the 
quantitative disagreement [45| between thermodynamic 
features of KCMs and real experiments is of little dy- 
namical consequence (4|| . In support of this perspective 
is the fact that KCMs have been remarkably successful 
in generating features of dynamical heterogeneity such 
as Stokes-Einstein decoupling, growing dynamical length 
scales, and excess tails in the real-space particle displace- 
ment function [H HH ES El • 

On the other hand, one may wonder if a deeper view- 
point would allow for an understanding of the kinetic 
rules that govern particle motion in the supercooled liq- 
uids. It is natural to speculate that such aspects might 
have roots in the thermodynamics of configurations. In- 
deed, simple local Monte Carlo "dynamics" can repro- 
duce all features of dynamical heterogeneity seen in New- 
tonian molecular dynamics simulations, and are based 
simply on making local moves that are configurationally 
allowed [4!|. Lattice models based on this concept are 
called "lattice glass models" (LGMs), and were first con- 
sidered by Biroli and Mezard [5(J. The rules for such 
models seem at first sight like that of KCMs. For ex- 
ample in the simplest versions of such models a particle 
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Figure 1: Comparison and distinction of a caricature of a 
kinetically constrained model with a lattice glass model. In 
the KCM any configuration is allowed, but move may only 
be made if a particle has at least one missing neighbor before 
and after the move. In the LGM, the global configuration is 
defined such that all particles must have at least one missing 
neighbor, and all dynamical moves must respect this rule. 
Note that the local environment around the moving particle 
is identical in this example, while the global configurations 
are distinct. Periodic boundary conditions are assumed for 
both panels. 



may move if it is surrounded by no more than a fixed 
number of nearest neighbors before and after the move 
[5ll I521 . 1531 HU . Locally this is identical to the type of dy- 
namical constraint that appears in the KCMs introduced 
by Kob and Andersen [43[ . However, this constraint must 
be met globally: all particles must have no more than a 
fixed number of nearest neighbors. As the density of the 
system increases, fewer and fewer configurations exist for 
which these constraints may be satisfied. It is thus the 
entropy of configurations that governs the slowing of dy- 
namics, intimately connecting the non-trivial thermody- 
namic weight of states accessible to the local dynamics. 
Indeed, LGMs can be solved exactly within the Bethe ap- 
proximation, or on Bethe lattices [j^JHH, and have been 
shown to have a glass transition due to the vanishing of 
the configurational entropy. The distinction between the 
KCM and LGM viewpoint is illustrated in Fig. [1] 

LGMs have been studied by a number of groups, but 
the focus has not generally been on real-space aspects 
of dynamical heterogeneity. For example, Coniglio and 
coworkers have developed a simple LGM that avoids crys- 
tallization and displays many features of typical glass- 
forming materials, including a growing multi-point sus- 
ceptibility (Xi(t)) [H,[54j]. On the other hand, this sys- 
tem appears to behave as a strong glass-former, with a 
stretching parameter close to one, and exhibits essen- 
tially no Stokes-Einstein violation. Our goal in this work 
is to survey in detail the dynamical behavior of a new 
LGM which bears similarity to the original Biroli-Mezard 
model but is not prone to crystallization. The main con- 
clusion that we draw is that LGMs are at least as realis- 
tic as KCMs in their description of all commonly studied 
features of dynamical heterogeneity. In this regard, sim- 
ple coarse-grained lattice models based on the thermody- 
namic weight of states are no less viable as fundamental 



caricatures of glassy liquids than are KCMs based on 
weights of trajectories. We conclude our work by high- 
lighting several key ways that LGMs and KCMs may be 
distinguished. We reserve the investigation of these com- 
parisons for a future study. Our paper is organized as 
follows: Sec. [II] outlines the model. Sec. IIIII discusses 
both simple averaged dynamics as well as aspects of dy- 
namical heterogeneity. In Sec. IIVI we conclude with a 
discussion of the meaning of our findings and the future 
directions to be pursued. 



II. MODEL 

Here we define the LGM that forms the basis of our 
simulations. The original model of Biroli-Mezard is quite 
prone to crystallization 50] . This fact makes its use prob- 
lematic for the study of glassy behavior since crystalliza- 
tion always intervenes before supercooling becomes sig- 
nificant. The crystallization problem persists on a square 
lattice for all binary mixtures we have studied. How- 
ever, we have found that certain generalizations of the 
Biroli-Mezard model with three species of particles are 
stable against crystallization for the densities that are 
sufficiently high that glassy dynamics may be clearly ob- 
served. 

Our model follows the original rules of the Biroli- 
Mezard model. Particles exist on a cubic periodic lat- 
tice of side L — 15 and each lattice site can contain only 
zero or one particle. All particles, at all times, must sat- 
isfy the condition a particle of type "m" must have m 
or fewer neighbors of any type. A neighbor is considered 
any particle in one of the 2d (^dimensionality) closest 
lattice sites along the cubic coordinate axes 1 . 

The particular three species model we employ is de- 
fined by 10% type 1 particles, 50% type 2 particles, 
and 40% type 3 particles. We denote this model the 
"tl54" model to indicate its basis in thermodynamics 
and to specify the types and percentages of each parti- 
cle. The composition of tl54 model was determined via 
trial and error by picking particle types with clashing 
crystallization motifs thereby frustrating crystallization. 
Crystallization was monitored by inspection of the an- 
gle resolved static structure factor, direct inspection of 



1 It should be noted that LGMs of the type described here in- 
volve extreme constraints that must be globally satisfied and 
are thus not realistic translations of off-lattice particle-based 
models. Such constraints might indeed induce artificial behav- 
ior, especially at higher densities. It would be most interesting 
to investigate "soft" versions of such models where constraints 
may be locally violated at the cost of an energy penalty. In 
this regard, such models would be the configurational analog of 
KCMs where dynamical constraints may be broken at the cost 
of an energy penalty, see Chandler, D. and J. P. Garrahan, "Dy- 
namics on the Way to Forming Glass: Bubbles in Space-time" 
arXiv:0908. 0418^1; Submitted to Annual Reviews of Physical 
Chemistry (2010). 
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configurations, and by monitoring bulk thermodynamic 
quantities. 

As discussed in the introduction, there appear to be 
strong similarities between the rules that govern KCMs 
such as the Kob-Andersen model and the tl54 model 
[43| . For example both models employ constraints with a 
maximum number of neighbors, but in the Kob-Andersen 
model this restriction only applies to the mobile parti- 
cles, while in the tl54 model applies to all particles. Our 
model does not require any special dynamics methods. 
We employ local canonical Monte Carlo "dynamics" via 
primitive translational moves (49J. Note that for the tl54 
model the energy can only be zero (no packing violations) 
or infinite (packing violation or overlap) , thus the accep- 
tance criteria reduces to rejection if there is a packing 
violation and acceptance otherwise. This allows us to 
implement an event-driven algorithm which accelerates 
the simulation of lattice dynamics [55| . 

For thermodynamic studies we employ grand-canonical 
Monte Carlo with both translational moves as well as 
particle insertion/deletion. Fig. [^contains a plot of the 
density of the system as a function of the chemical po- 
tential of type 1 particles. Models which crystallize (such 
as original binary model of Biroli-Mezard) have a sharp 
jump in this curve at the crystallization point. Clearly, 
this feature is absent in the tl54 model. For comparison, 
both curves are displayed 2 . 



III. DYNAMICAL BEHAVIOR 
A. Simple Bulk Dynamics 
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Figure 2: Crystallization thermodynamics in LGM. Top: The 
tl54 model, ^ti refers to the chemical potential of the type 
1 particles. The maximum density observed for the 15 3 lat- 
tice is .5479 (exactly 1849 out of 3375 lattice sites occupied). 
The three plotted quenching rates vary between a .01 and .05 
increase of /xi per 10000 cycles. Bottom: A close up of the 
equivalent plot for the BM model. Note the clear disconti- 
nuity upon crystallization. Slower /x-increase rates produce a 
sharper discontinuity. 



In this subsection we describe the behavior of a simple 
2-point observable, namely the self-intermediate scatter- 
ing function [58|, defined as 

F s (k,t) = (±Y,^ lrdt) - r * m )- (1) 



2 A subtle issue arises in the nature of glassy behavior observed 
in the tl54 model outlined in this work. LGMs could have a 
dynamical percolation-like transition, as in the spiral model [56l . 
This has been indeed found in some LGMs on the Bethe lattice 
[57j and would slow down the dynamics for reasons completely 
different from the diminishing of the configurational entropy. If 
there is a low-lying crystal phase then one can show that this 
dynamical percolation-like transition cannot take place in finite 
dimension. Although we have not found a crystal phase for the 
model, the existence of such a transition seems unlikely and irrel- 
evant for our present work. First, it can be shown that blocked 
structures, if they exist, have to verify much more constraints 
than in the spiral model [5611 . Second, we have found that the re- 
laxation time growth of the persistence functions with increasing 
density in local canonical Monte Carlo simulations are similar 
to those under grand-canonical dynamics, which cannot contain 
any blocked structure. The union of these two facts render the 
dynamical blocking scenario highly unlikely. 



We measure F s (k,t) only for the type-2 particles which 
are present in the greatest fraction for the three distinct 
species. Throughout this paper, we report fc-vectors us- 
ing k' , where k — ^-k' '. We have checked that F s (k,t) 
is qualitatively similar for the other species of particles. 
The relaxation of F s (k, t) of the system at the wavevector 
k' = 5 (k = ^p) for various densities is shown in Fig. [3] 
The bulk of the decay may be fit to a stretched exponen- 
tial function, F s (k,t) = exp(-(t/r Q (fc)) /3 ( fc ^). As is cus- 
tomary, the alpha-relaxation time is found by the value 
F a (T a ) — 1/e and the f3{k) exponent is determined by a 
direct fit to the terminal decay. We find that for densities 
below approximately p — 0.48 the value of (3 saturates at 
the expected value (3 = 1 characteristic of simple non- 
glassy dynamics, while for the highest density simulated, 
[3 = 0.7. This behavior, over a similar range of super- 
cooling, is reminiscent of the behavior found in atomistic 
models of glass- forming liquids [59t l6fj • In order to bet- 
ter reveal the relaxation behavior, F s (t) is also displayed 
on a log-log vs. log-time scale. In this plot, the slope 
of the long time growth is related to the exponent (3. We 
have found that the values of (3 extracted from the slopes 
of the long time portion of the log- log vs. log plot indeed 
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coincide with that found by a direct fit to a stretched 
exponential form. At the highest densities a shoulder 
appears in the short time relaxation. This feature is in- 
dicative of a secondary relaxation feature perhaps akin 
to beta-relaxation in realistic glass- forming liquids. It 
should be noted, however, that the amplitude of this fea- 
ture is very close to unity. This is quantitatively distinct 
from the plateau values expected in atomistic off-lattice 
models [53, [6(| and even LGMs with more complicated 
lattice degrees of freedom [H, H3| , but is similar to that 
encountered in simple spin models such as variants of the 
Random Orthogonal Model (6lT |. 

As is typical of fragile glass-forming systems, the tl54 
model exhibits relaxation times that do not follow the 
(generalized) Arrhenius form [62] . This behavior is illus- 
trated in Fig. |U At low densities, plots of log(r) versus 
p indeed follow a straight line, however in the vicinity of 
p ~ 0.5 the plot of r versus p deviates from this straight 
line and the functional density dependence of the relax- 
ation time becomes much more precipitous. While we 
have not attempted to quantitatively characterize this 
density dependence, it should be noted that the onset of 
increased sensitivity to changes in density occurs is the 
same narrow window that marks the noticeable decrease 
in the values of the stretching exponent (3. 



B. Motion on the Atomic Scale 

We begin our discussion of the nature of heterogeneous 
dynamical behavior in the tl54 LGM by observing the 
qualitative details of particle motion under supercooled 
conditions. This will set the stage for analysis of quanti- 
tative measures of dynamical heterogeneity in the model. 
For the sake of comparison, we also investigate the analo- 
gous behavior in the Kob- Andersen model. This compar- 
ison is useful because it suggests how models with similar 
local rules but different global rules (rooted in either the 
purely kinetic or thermodynamic basis of the particular 
model) may give rise to distinct dynamics at the particle 
scale. 

We start by simply observing the patterns of mobil- 
ity in real space starting from a set initial condition of 
the tl54 model found at a given density after equilibra- 
tion. A similar analysis has been performed recently by 
Chaudhuri et al. for the Kob- Andersen model, where no 
equilibration is required since all initial configurations 
with a set density of defects are allowed [63| . For a theo- 
retical description of the dynamics of the Kob- Andersen 
model, see [64j. We note that, as expected, the tl54 
model exhibits regions of spatially localized particle ac- 
tivity against a backdrop of transiently immobilized par- 
ticles. A rather remarkable feature of the patterns of 
mobility in this model is that we find evidence of string- 
like motion, where a group of particles moves over a short 
distance, each taking the place of the previous particle in 
the string [(55|, [66|. This motif can be seen mostly on 
timescales less than the a-relaxation time, but occasion- 
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Figure 3: Decay of the self-intermediate scattering function 
F s (k,t) for k' = 5 (fc = ^-k'). Densities are .3, .4, .45, .48, .50, 
.51, .52, .53, .535, .5375, .5400, .5425 from fastest relaxation to 
slowest relaxation. These densities are used in all plots in this 
paper unless otherwise indicated. Top: Plotted on a linear- 
log scale. Bottom: Same data as upper panel plotted on a 
log(— \og 10 (F s (k, t))) vs log(t) scale. Lowest density curves 
are at the top left. 



ally string-like motion may be seen to persist on longer 
timescales. This behavior is demonstrated in Fig. 

The behavior of particle motion observed in the Kob- 
Andersen model is somewhat different than that seen 
in the tl54 model as described above. As in the tl54 
model, and as observed by Chaudhuri et al., motion in 
the Kob- Andersen model shows similar activity regions 
in the vicinity of defect sites giving rise to heterogeneous 
motion. However, the boundaries between active and 
inactive regions at comparable timescales appear to be 
more distinct in the Kob- Andersen model. Furthermore, 
the particle scale motion in the Kob- Andersen is much 
more isotropic, exhibiting much fewer cases of directional 
mobility compared with the tl54 model. It would be in- 
teresting to compare the two models by quantifying this 
difference via the type of directional multi-point correla- 
tors devised by Doliwa and Heuer [U- It is not clear if 
the difference between the models is related to the funda- 
mental distinction between LGMs and KCMs or just the 
specifics of the particular models considered. In particu- 
lar, the tl54 is a multi-component model, unlike the Kob- 
Andersen model. The string-like motion on short time 
scales seems to occur predominantly on the rather rough 
boundaries of slow clusters (68j. This behavior, reminis- 
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Figure 4: Top: r a (time at which F a (k,t) = 1/e) as a func- 
tion of density, p. Plotted for k' = 1, 2, 3, 4, 5, 6, 7, with lowest 
k at the top. Center: Beta stretching exponent of F 3 (k,t) 
(from terminal fits F 3 (k,t) ~ exp( — (t/T a )' 3 )). Lowest curve 
is at the top of the plot. Bottom: Plot of log scale r a against 
chemical potential p of type 2 particles. The behavior is con- 
sistent with Tc = 5.7exp(— 21/i2/(/^2 — 24)). 




Figure 5: Examples of string-like motion apparent in the tl54 
model, (a) An example of a string with all neighboring par- 
ticles removed, (b) A similar string in the context of other 
particles. Note that here the string is truly isolated in space, 
away from other mobile particles. In these figures, type 1 
particles are white, type 2 particles are blue, and type 3 are 
green. Sites occupied at the initial time but vacated at the 
final time are shown in red. These pictures show only the 
differences in position of particles between the origin of time 
and the final time, not the path the particles took to achieve 
that displacement. All figures are at a density of .5400, with 
At times in (a) 251, (b) 199526. The a-relaxation time for 
k' = 5 at this density is about 7.8 x 10 6 




(a) (b) 

Figure 6: Examples cluster shapes in the (a) the tl54, model, 
density p — .5400 and (b) the Kob-Andersen model, density 
p = .8500. Arrows indicate motion between initial and final 
times. Time separation is l/10th of the a-relaxation time. In 
the tl54 model, we see more fractal and disconnected clusters, 
while in the KA model, mobile domains tend to be smoother 
clusters. 



cent of the picture of dynamic heterogeneity that put 
forward by Stillinger [691 ] . might be strongly influenced 
by compositional heterogeneity. A useful way to address 
general issues related to how the initial configuration con- 
strains subsequent dynamics would be a systematic iso- 
configurational ensemble analysis comparing LGMs and 
KCMs [70]. This will be the topic of a future publication 

a 

In the next few sections, we discuss how some of the 
most important indicators of dynamical heterogeneity in 
supercooled liquids manifest in the tl54 model. The 
quantities that we discuss are the magnitude of violations 
of the Stokes-Einstein relation, exponential tails (indica- 



tive of hopping transport) in the van Hove function, the 
existence of a Fickian length scale and the development 
of a dynamical length scale quantified by the multi-point 
function S4(q,t). Unless otherwise stated, specific corre- 
lation functions and transport coefficients are calculated 
with respect to type-2 particles. 



C. Stokes-Einstein Violation 

In typical fluids a mean-field linear-response relation- 
ship asserts that the product of the tracer particle dif- 
fusion constant and the fluid viscosity divided by the 
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temperature is a constant 58]. This connection be- 
tween diffusion and dissipation is known as the Stokes- 
Einstein relationship, and empirically is known to hold 
even at the atomic scale in liquids over a wide range 
of densities and temperatures. In supercooled liquids, 
the Stokes-Einstein relation generally does not hold 
[U HH H3, [H, [M HE M, UM- In fact, the product 
of the diffusion constant and the viscosity of a liquid 
may exceed that expected from the Stokes-Einstein re- 
lation by several orders of magnitude close to the glass 
transition. There are many theoretical explanations for 
Stokes-Einstein violations in supercooled liquids, which 
essentially all invoke dynamical heterogeneity as the fun- 
damental factor leading to the breakdown of the simple 
relationship between diffusion and viscosity. It should be 
noted that similar relationships hold between the diffu- 
sion constant and the self and collective time constants 
associated with the decay of density fluctuations. In 
this work we focus on the relaxation time of the self- 
intermediate scattering function defined above as our 
proxy for the fluid viscosity. 

It is well known that the product Dr a , where r Q is 
the a-relaxation time of the self-intermediate scattering 
function shows a strong temperature/density dependence 
in both realistic atomic simulations as well as in the class 
of KCMs that describe fragile glass-forming liquids. No 
direct studies of this quantity have been made in LGMs. 
The LGM of Coniglio and coworkers would appear to 
show essentially no Stokes-Einstein violations because 
the diffusion constant and the relaxation time may both 
be fit to power laws with exponents that have, within 
numerical accuracy, the same magnitude 1541]. This, 
however is not surprising since many of the features of 
the model resemble those of a strong glass-forming sys- 
tem, where violations of the Stokes-Einstein relation are, 
at most, weak. The features of the tl54 model with re- 
gard to non-exponential relaxation and the density de- 
pendence of the relaxation time r Q indicate that this 
model behaves more like a fragile glass former. Thus, 
we expect clear violations of the Stokes-Einstein relation. 
Indeed, as shown in Fig. [7J Dr a increases markedly as 
density is increased. Over the range densities that we can 
access, the magnitude of the violation is very similar to 
that seen in the canonical Kob- Andersen Lennard Jones 
mixture over a comparable range of changes in relax- 
ation time [13] • Interestingly, violations begin to become 
pronounced at densities similar to where the relaxation 
times and stretching exponents become strongly sensitive 
to increased density. Thus, a consistent onset density is 
observed as in more realistic atomistic systems. 



D. van Hove Function 

It is now rather well established that an additional 
"quasi-universal" feature of dynamical heterogeneity 
near the glass transition is contained in the shape of the 
real-space van Hove function [1|| H3, [H, Hsj] . In particu- 




0.25 0.3 0.35 0.4 0.45 0.5 



p (total) 

Figure 7: Violation of the Stokes-Einstein relation, Dr a ~ 
constant, using r a at k'—5. Data has been normalized to 
Dr a — 1 at the lowest density. 

lar it has been argued the tails of the self van Hove func- 
tion should be approximately exponential in form. These 
"fat tails" imply that the rare particles that do undergo 
large displacements exist in populations in excess of what 
would be expected in a purely Gaussian displacement dis- 
tribution. While non-Gaussian tails should be expected 
of any distribution for the wings that fall outside of limits 
of bounds set by the Central Limit Theorem, the palpable 
exponential tails in supercooled liquids imply large non- 
Gaussian effects indicative of transport that is strongly 
effected by heterogeneous hopping motion. 

Here, we demonstrate that such effects occur in the 
tl54 model in a manner similar to that seen both in ex- 
periments in colloidal and granular systems as well as in 
computer simulations of atomic systems. Fig. [8] shows 
the self part of the real-space van Hove function, 

G s (x,t) = (S(x-\Z-(r i (t)-r i (0))\)), (2) 

for the type two particles in the tl54 model. Because we 
are on a lattice, we restrict our distances along the three 
coordinate axes x individually in our calculation. We 
see that for times of the order of the a-relaxation time, 
these tails are clearly visible. For very long or short time 
scales, the shape of the tail deviates somewhat from the 
more exponential form exhibited at intermediate times. 
This behavior is quite similar to that seen in simulations 
of atomistic systems [13, H|| , and is fully consistent with 
the behavior found in KCMs [Hj]. 

E. Fickian Length 

Related to the existence of excess tails in the van Hove 
function is the existence of a length scale that character- 
izes the anomalous transport. More specifically, the ex- 
ponential tails in the van Hove function are distinguished 
from the Gaussian form of the displacement distribution 
obtained at relatively short distances for fixed times. The 
crossover from Fickian to non-Fickian behavior should be 
characterized by time scales as well as length scales over 
which this crossover occurs. A non-Fickian length scale 
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Figure 8: van Hove function for p = .5375 and various times. 
Distances are measured independently along each coordinate 
axis. The times plotted, from left to right, are 10 5 , 316227 
(approx. the a-relaxation time), and 10 6 . An exponential fit 
to the tail of the t = 316227 case is shown by a dotted line. 
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Figure 9: fc-dependent diffusion constant = 1/ (fc' 2 T a (fc)). 
Densities of .3000 (upper) and .5425 (lower). The higher den- 
sity curve is multiplied by a scale factor of 2.992 x 10 5 for ease 
of comparison. A dotted flat line is included for reference of 
behavior expected in the purely Fickian case. 



may be defined by examining the fc-dependent diffusion 
constant D(k) = -^-^ [30l. I3U |32| . The wavevector that 
characterizes the crossover from the expected diffusive 
behavior to an anomalous regime is inversely related to 
such a length scale. In Fig. [5] we plot D(k'). Clearly, 
as the density is increased, the length scale separating 
the Fickian and non-Fickian regimes increases. This be- 
havior is consistent with that found in KCMs and sim- 
ulations of atomistic glass-forming liquids. It should be 
noted that Stokes-Einstein violations, the development 
of exponential tails in the self van Hove function, and a 
well-developed Fickian length scale are all manifestations 
of related aspects of dynamically heterogeneous motion 
in supercooled liquids [271 ] . 



F. x-i and Si Fluctuation Measures 

The Fickian length scale is merely one length scale that 
arises naturally in systems where dynamics become in- 
creasingly heterogeneous. Perhaps more fundamental is 
the growth of dynamical length scales associated with 
multi-point correlations of the dynamics. Supercooled 



liquids do not show simple static correlations that would 
indicate a growing correlation length. It should be noted 
that this does not exclude growing static correlations of a 
more c omp lex kind, for example point-to-set correlations 
[3, [H, [74|. Regardless, cooperativity in dynamics may 
be measured via first defining a local overlap function 

= ^E^""^ [«»(k. (Ari(t))) - F s (k,t)} 

i 

(3) 

where Ar^(t) = r^O) — r^t). fk(q,t) is defined for one 
configuration, and the average is over all k and q with 
the magnitudes k and q. Then, S^q) is defined as 

S i (q) = N(\6f k (q,t)\ 2 ) (4) 

where this average is over the most general ensemble of 
configurations [361 ]. The X4 value is defined as the limit 
S4(q — > 0). Xi(t) ma Y be calculated strictly at q = 
from 

Xi(t)=N(j6f k (q = 0,t)\ 2 ) (5) 

where the average is over the entire ensemble and all 
k consistent with the magnitude of k. Note that, as dis- 
cussed in (36j , the value of Xi{t) computed in this manner 
is a lower bound for the extrapolation of Si(q — > 0,t). 

The quantity S4,(q,t) is a multi-point dynamical ana- 
log of S(q). Just as the low q behavior of S(q) indicates 
a growing (static) length scale in systems approaching a 
second order phase transition, scattering from dynami- 
cally heterogeneous regions undergoing cooperative mo- 
tion will manifest growth in the amplitude of the low q 
region of Sf (q,t), indicative of the size scale of the dy- 
namical correlations for systems approaching the glass 
transition. 

The behavior of the quantity 5| (q, t) is shown in Fig. 
[TU1 Only type-2 particles have been used in the calcula- 
tion. As can clearly be seen, for densities above p ~ 0.5 
which constitutes the onset density of this system, the 
low q behavior shows a marked upturn as q — > 0. The 
growth of S4 (q,t) as q — > as density is increased sug- 
gests a growing length scale as supercooling progresses. 
This non-trivial behavior is what is found in atomistic 
simulated systems. Future work will be devoted to a pre- 
cise characterization of the length scale that may be ex- 
tracted from £4' (q, t) in the tl54 model so that a compar- 
ison may be made with recent work detailing the behavior 
of this length in realistic off-lattice systems [H, [75[ . 

IV. CONCLUSION 

In this paper we have presented a new LGM based on 
the original Biroli-Mezard model 50] . Via the introduc- 
tion of an additional species of particle, we have demon- 
strated that our model is stable against crystallization. 
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Figure 10: Top: Top:Plot of Si(q,t) at r a for densities 0.51, 
0.52, 0.53 and 0.54. Bottom: Plot of Xi(t) for the same 
densities. Peak values correspond to lower bounds for of the 
value of S4,(q,t) in the upper panel at q = 0. 

This fact allows us to study sufficiently high density con- 
figurations that manifest features of dynamical hetero- 
geneity. Unlike some previous LGMs, our model exhibits 
the canonical features of a fragile glass-former. In terms 
of the gross features of relaxation behavior, our LGM 
shows behavior similar to the standard Kob-Andersen 
Lennard-Jones (KALJ) mixture. In particular, we find 
that the degree of violation of the Stokes-Einstein rela- 
tion and the magnitude of stretching in the decay of the 
self-intermediate scattering function track the relaxation 
times at densities above the onset of supercooling in a 
manner consistent with that seen in the KALJ system. 
Features of dynamical heterogeneity such as exponential 
tails in the van Hove function, the growth of a dynamical 
length scale as quantified by the function Si(q, t), Stokes- 
Einstein violations and the emergence of a Fickian length 
scale all occur in a manner expected from experiments 
and simulations of fragile glass-forming liquids. 

The similarity between the description of dynamic het- 
erogeneity found in KCMs and LGMs stands in stark con- 
trast to the underlying foundations of the models them- 
selves. As emphasized in the introduction, KCMs are 
based on a constrained dynamics for which the num- 
ber of available dynamical paths leading to relaxation 
becomes increasingly rare as the density increases and 
the number of defects decrease. In KCMs all real-space 
configurations at a fixed number of defects (excluding 
rare blocked configurations) are equally likely. On the 
other hand LGMs are based on transitions between real- 



space configurations that become increasingly scarce as 
the density is increased. This is not to say that there 
is not a facilitated-like dynamics in LGMs. On the con- 
trary, as we have demonstrated in sec. IIII1 local and 
sometimes anisotropic dynamics may be generated nat- 
urally in LGMs without the explicit introduction of fa- 
cilitating defects. An important message that emerges 
from this study is that the phenomenology of dynamic 
heterogeneity is not sufficient to distinguish pictures or 
validate models based on transitions between sets of states 
in configuration space from those based on sets of paths 
in space-time. 

How then might these pictures be differentiated? 
While contrasting competing models that generate seem- 
ingly similar dynamical behavior is a difficult endeavor, 
several possible studies might be useful for this task. Here 
we outline four avenues that could provide key informa- 
tion that distinguish the purely dynamical picture from 
one based on transitions thermodynamic states. 

a) The mosaic length scale: The Random First Order 
Theory (RFOT) of Wolynes and coworkers posits the ex- 
istence of a static length scale which is defined by the 
region over which particles are pinned by the surround- 
ing self-generated amorphous configuration [H, [12, l38| - 
This length scale also exists in KCMs, but it is decou- 
pled from the relaxation dynamics of the system [7t| . 
Recent atomistic computer simulations have successfully 
located the mosaic length scale [38]. It would be quite 
useful to perform an analysis similar to that devised by 
Jack and Garrahan for LGMs [76] . Since LGMs are based 
on the entropy of real-space configurations, it is expected 
that here the mosaic length does couple to the glassy dy- 
namics. Since LGMs are much simpler than atomistic 
off-lattice models, the direct study of the mosaic length 
(and point-to-set correlations in general) in LGMs might 
provide key avenues for the testing of the putative cou- 
pling between relaxation and such length scales in simu- 
lated atomistic systems. 

b) Correlations between configurational entropy and 
dynamics: Empirical correlations between the configu- 
rational entropy and the a-relaxation time have been 
noted for many years, and this correlation lies at the 
heart of several prominent theories. Such correlations are 
still widely debated, but seem to hold at least crudely in 
many glass-forming systems [z3, [78[. LGMs should be 
expected to exhibit such correlations, while it is known 
that KCMs do not exhibit such correlations. Recently 
Karmakar et al. purported to show that finite-size effects 
of the a relaxation time follow precisely the Adam-Gibbs 
relation between the configurational entropy and the en- 
relaxation time in the KALJ system [39(. If true, such 
correlations would be a challenge to KCMs, since it is dif- 
ficult to envision how the configurational entropy would 
track the a-relaxation time for different system sizes if it 
were not a crucial component of relaxation phenomena. 
Such correlations, however, are subtle to measure since 
the Adam-Gibbs relationship is an exponential one and 
the apparent correlation could depend on the somewhat 
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indirect computational method used in [39j to define the 
configurational entropy. It would be most useful to in- 
vestigate such effects in the simpler LGMs, which might 
provide a cleaner means of isolating the configurational 
entropy. It should be noted that finite size effects do 
appear to follow an approximate Adam-Gibbs relation- 
ship in at least one other lattice model [7^|. Such stud- 
ies might spur more detailed investigations in simulated 
atomistic systems thus allowing for a clear comparison 
between LGMs, KCMs and more realistic systems. 

c) Single-particle and collective predictability ratios: In 
an important piece of work, Jack and Berthier devised 
metrics that access the degree to which single particle 
and collective dynamics are deterministically predicted 
by a set initial configuration over a given time scale [80l | . 
KCMs and LGMs differ in how allowed configurations are 
constructed. KCMs have explicit defects, while configu- 
rations in LGMs are determined by global constraints, 
and thus do not contain explicit defects. Since the very 
composition of initial conditions differ markedly in these 
models, one expects that the metrics defined by Jack and 
Berthier would behave differently in KCMs and LGMs. 
Thus, it would be very profitable to examine the density 
and temperature dependence of the single particle and 
collective predictability ratios in KCMs and LGMs as 
a possible means of distinguishing between state-based, 
and dynamical constraint-based pictures [7l| . 

d) Evolution of the facilitation mechanism approaching 
the glass transition: Although in both KCM and LGM 
pictures facilitation plays an important role in the relax- 



ation of the system, a peculiar and different temperature 
and density evolution is expected. In particular, in the 
KCM picture, facilitation is due to the motion of mobility 
regions or defects. Dynamics slows down, and concomi- 
tantly dynamic heterogeneity increases, because these re- 
gions become rarer approaching the glass transition. A 
crucial assumption is that these defects are conserved or 
at least that non-conservation is a rare event that be- 
comes rarer at lower temperature/high density. These 
assumptions impose important constraints on the evolu- 
tion of the facilitation mechanism. Thus, it would be 
very interesting to examine this issue for example using 
the cluster analysis developed in [sij] to study the relax- 
ation dynamics of granular systems. 

Investigation of these and other studies aimed at 
distinguishing the underlying pictures that LGMs and 
KCMs are based on will be the subject of future work. 
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